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Abstract 

The existence of bulk sterile neutrinos in theories with large extra dimen- 
sions can naturally explain small 4-dimensional Dirac masses for the active 
neutrinos. We study a model with 3 bulk neutrinos and derive various con- 
straints on the size of the extra dimensions from neutrino oscillation experi- 
ments. Our analysis includes recent solar and atmospheric data from SNO and 
Super-Kamiokande, respectively, as well as various reactor and accelerator ex- 
periments. In addition, we comment on possible extensions of the model that 
could accommodate the LSND results, using natural parameters only. Our 
most conservative bound, obtained from the atmospheric data in the hierar- 
chical mass scheme, constrains the largest extra dimension to have a radius 
R < 0.82 ^m. Thus, in the context of the model studied here, future gravita- 
tional experiments are unlikely to observe the effects of extra dimensions. 



*Pcrmancnt address: Department of Physics and Astronomy, University of Pennsylvania, 
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1 Introduction 



The substantial body of theoretical work done on models with Large Extra Dimensions 
(LED's) has provided us with new insights on some old problems. Even though the 
original motivations were mostly related to the question of gauge hierarchy , various 
other applications have been found and studied. Among these is an explanation 
for small but non- vanishing neutrino masses 0. In the more traditional scenario 
with a large energy desert between the scales of electroweak symmetry breaking and 
grand unification, the small Majorana masses can be easily generated through the 
see-saw mechanism. In the models with LED's, however, this is not possible, since 
physics at energy scales well above a TeV is no longer described by a renormalizable 
quantum field theory. There is, however, an alternative mechanism. In these models, 
all the Standard Model (SM) particles, including left-handed neutrinos, have to be 
confined to a 4-dimensional subspace (3-brane) inside the full space-time. On the 
other hand, if an SM singlet fermion (such as the right-handed neutrino) is present, it 
can propagate in more than four dimensions. The large volume of the extra dimensions 
leads to a suppression of the wave function of this fermion on the brane, which in turn 
allows small Dirac neutrino masses to be generated naturally. This point contributes 
positively towards a comparison of the models with LED's with those of the traditional 
energy desert paradigm. 

From the four- dimensional point of view, a higher- dimensional SM singlet fermion 
can be decomposed into a tower of Kaluza-Klein (KK) excitations. These states do not 
have SM charges, and can therefore be classified as "sterile" . However, the KK states 
do not completely decouple from the system: there are mixings between them and 
the lowest-lying, active neutrinos. Thus, they can participate in neutrino oscillations, 
acting effectively as a large number of sterile neutrinos. The implications of this 
picture have been studied in a number of papers 0, H, ^, ^, |[|]. 

The recent experimental results of the SNO collaboration |^, in conjunction with 
the data from Super-Kamiokande , have yielded strong constraints on the contribu- 
tion of any sterile state to the solar neutrino anomaly. Motivated by the new data, we 
have reexamined the phenomenological constraints on models with extra-dimensional 

t There have been a number of other models which combined the ideas of LED with additional 
ingredients, such as small Majorana masses for the brane neutrinos 0. 
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neutrinos. We have considered the bounds coming from reactor and accelerator ex- 
periments, as well as measurements of fluxes of solar and atmospheric neutrinos. We 
have studied a model 0, ^ with 3 bulk neutrinos which have Yukawa couplings to 
the 3 active brane neutrinos. We will refer to this model as the (3, 3) model, where 
the first integer corresponds to the number of brane neutrinos and the second is 
the number of bulk neutrinos. Unlike most previous studies, we are considering the 
situation in which the dominant effects for the solar and atmospheric neutrinos are 
oscillations amongst the three active states, with oscillations into the sterile KK ex- 
citations regarded as a perturbation to be constrained by the data. This framework, 
which is motivated by the SNO and Super-Kamiokande data, allows us to obtain 
simple analytic results throughout the analysis. We emphasize that the principal role 
of the extra dimensions in this case is to provide a natural framework for small Dirac 
neutrino masses. 

In this paper, we will assume that the manifold on which the extra dimensions 
are compactified is highly asymmetric, with one dimension much larger than the rest. 
Unlike most collider, astrophysical and cosmological bounds, which are largely in- 
dependent of the compactification manifold's shape and only constrain its volume, 
the analysis of this paper will yield bounds on the size of the largest dimension. The 
bounds we obtain are sufficiently tight to rule out (in the context of the model consid- 
ered here) effectively the possibility that the near-future Cavendish type experiments 
will observe a signal for large extra dimensions. Although this conclusion clearly de- 
pends on the assumed properties of the neutrino sector of the theory, we note that 
the (3, 3) model seems to be the simplest and most natural way to incorporate small 
neutrino masses in theories with LED's. 

Along with the solar and atmospheric neutrino anomalies, evidence for neutrino 



oscillations has been reported by the LSND experiment |10[. The LSND results cannot 
be accommodated within the minimal (3, 3) model without the introduction of ad 
hoc parameters, such as brane Majorana masses. We propose two extensions of the 
(3, 3) scheme which can address these results, using only the natural parameters of 
the model. The first proposal introduces two new sterile states, one on the brane and 
the other in the bulk. Thus, this model is of the (4, 4) type. The second proposal uses 
an additional LED of a different radius, in which a sterile bulk state is sequestered. 
The paper is organized as follows. In the next section, we describe our setup and 
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introduce the necessary notation. The motivation behind our analysis as well as our 
approach and choice of parameters are discussed in section |^. Section | discusses 
constraints on the model from reactor, accelerator, atmospheric, and solar data. We 
present our results in the context of the hierarchical, inverted, and degenerate mass 
schemes. In section ^, we address the LSND results and outline the aforementioned 
proposals for accommodating them. Our conclusions are presented in section ^. 

2 Formalism and Notation 
2.1 General Setup 

As required in theories with LED's, we will assume that the Standard Model (SM) 
fields, including the three families of left-handed neutrinos z/£ {a = e, fi, r) and the 
Higgs doublet H, are confined to a four- dimensional brane, while gravitational fields 
propagate in a space-time with 6 >2 additional compactified dimensions of volume Vs. 
The fundamental gravitational scale of the higher-dimensional theory Mp is generally 
taken to be close to the weak scale to solve the hierarchy problem. To reproduce the 
measured strength of four-dimensional gravity at large distances, we require 

Ml, = Mp-^Vs, (1) 

where Mpi ^ 2 x 10^*^ GeV is the reduced Planck mass. 

We postulate the existence of three families of SM singlet bulk fermions which 
can propagate in 4-1-5 dimensions. Yukawa couplings of \E'" to the left-handed neutri- 
nos on the brane give rise to masses for active neutrino species which are naturally 
of the right size to accommodate the solar and atmospheric neutrino anomalies. The 
same couplings lead to mixings between the active species and the higher Kaluza- 
Klein (KK) components of the bulk fermions, which from the four-dimensional point 
of view appear as sterile neutrinos. 

Throughout this paper, we will assume that one of the dimensions is compactified 
on a circle of radius R which is much larger than the sizes of the other dimensions. In 
this case, only the KK excitations of the bulk neutrinos corresponding to the largest 
dimension will be relevant at low energies, and our treatment will be essentially five- 
dimensional. Apart from its simplicity, this model is interesting for the following 
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reason. LED's could be discovered by Cavendish type experiments which search 
for deviations of the gravitational force from Newton's law at short distances. At 
present, these experiments probe distances of order 0.2 mm |Tl|, and sensitivities 



of order 0.05 mm can be achieved in the near future [0]. If the compactification 
manifold is symmetric, astrophysical constraints on the radii of the extra dimensions 



imply that these sensitivities will not be sufficient to detect a signal |T3[. However, 
the astrophysical and other high-energy constraints primarily restrict the volume of 
the extra dimensions, while the Cavendish type experiments are sensitive to the size 
of the largest dimension. Thus, in models with highly asymmetric compactifications, 
such as the one we are considering, these experiments still have a chance of observing 
deviations from Newton's law. Unfortunately, the constraints obtained in this paper 
indicate that even this possibility may be already ruled out, if the neutrino sector of 
the model has the properties assumed in our study. 

As we will see below, the simple setup described here cannot accommodate the 
positive result reported by the LSND neutrino oscillation experiment. In section 5, we 
will propose two simple extensions of this setup that might be capable of explaining 
the LSND result. 



2.2 The (3, 3) Model 

From the 4-dimensional point of view, the 5-dimensional fermion \I' can be decom- 
posed into two Weyl fermions, ipL and ipR- The action of the model is given by 

S = J d^xdy i'^'^FA^^^" + J d^x (tV'l^f.d^'ul + XapHvl^'^{x,0) + h.c.) , (2) 

where F^, A = 0, . . . ,A are the 5-dimensional Dirac matrices. Note that this action 
conserves lepton number; in particular, we do not introduce Majorana masses for 
the left-handed neutrinos. The Yukawa couplings Xa/3 have dimensions of (mass)"'^''^. 
Since the only mass scale in the problem at this point is Mp, we introduce dimen- 
sionless Yukawa couplings via 

Kp = KpM'J\ (3) 

We will assume that hap are of order one. 

Let us decompose the 5-dimensional fermions '4'l,r into a tower of KK states 
V'i"]?; = . . . oo. It turns out that certain linear combinations of the KK states 
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are completely decoupled from the left-handed neutrinos, and therefore need not be 
considered. The states that do couple are given by 



«(0) 






a(n) 
^R 


= n=l...oo; 






= ^(^?"^+^?-"^), n=l...oo. 


(4) 



In this notation, the mass terms resulting from take the form 

(oo \ oo 

17^ + 72 + E + h.C. (5) 

n=l / n=l ^ 

The Dirac mass matrix is given by 

= h^pivMp/Mpi), (6) 
where v = {H) = 246 GeV, and we have used Eq. (g), (MIKs)!/^ ^ Mp^/Mp 

-I /ty 

(note that the KK modes have a prefactor proportional to V"^ ), and Eqs. If 
Mp/Mpi <^ 1, as is natural in models with LED's, the common scale of Dirac masses 
is well below the electroweak symmetry breaking scale v even for order-one Yukawa 
couplings. For example, with Mp ~ 100 TeV, and 0.1 ^ hajs ^ 1, Dirac masses of 
the active species will coincide with those required by the solar and atmospheric data. 
This is the higher-dimensional version of the see-saw mechanism. 

Diagonalization of the mass terms in can be performed in two steps. First, we 
find 3x3 matrices / and r which diagonalize the Dirac mass matrix, = r'^m^l = 
diag (mf, mf*, mf). We define 

















= 0. , 


. . oo; 








n = 1 . . 


. . oo. 


(7) 



The transformation of the ul by r"* greatly simplifies the problem by leaving the 
second term in (1) diagonal. Then, (|]) takes the form 

j:v'^M,u',^ + h.c. (8) 

i=l 
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where z/£* 



,/i ,/«(2) ^ 

. l^j. 1 • • - J 



I v'^^^^ , z/^'*'^'' , v'^^'^^ 1 ■ ■ - Y" 1 and Mj is an infinite- 



dimensional matrix given by 



m, 



D 



D 



V2m- 





1/R 







2/R 



V 



(9) 



To complete the diagonalization, we need to find infinite-dimensional matrices L 
and R such that R^MiLi is diagonal; then, the mass eigenstates are given by z/)^ = 
L|z/^\ z>)j = R\v'^. In this paper, we are primarily interested in the decomposition of 
the fiavor-basis brane neutrino states z/£ in terms of the mass eigenstates. These are 
given by 

On 



yl = Y^rY.L>r- (10) 

i=l n=0 

The easiest way to find L is to observe that it has to diagonalize the Hermitian matrix 
M^M, and therefore consists of its eigenvectors. This procedure was performed in 
Refs. 0, H, H. The component of L that enters Eq. (p!0D is given by 

2 



0n\2 
i J 



1 + 7r^a/2 + 2\t^'/i, 



(11) 



where = y/2mfR, and Aj""*^ are the eigenvalues of the matrix R^M^M, which are 
found from the equation 



A,, 



TT 



cot(7rAi 



0. 



The mass of the mode z/]^^"'' is simply given by Xf^' /R. Similarly, 

jkn _ T On 



An) 



(12) 



(13) 



where = 1 . . . oo and n = . . . oo. 

Since our philosophy will be to consider the effects of extra dimensions as small 
perturbations on top of the oscillations amongst zero-modes, we will be mostly inter- 
ested in the case where R~^ ^ mf , that is, the limit of small ^j. In this limit, we 
find 



Af^ = k + ^ij + ..., = | + ... 

* 2A; ' * A; 

and lI'^ = 0{^f), where j ^ k = 1 . . . oo. A quantity of crucial interest is the 
probabihty of finding a neutrino of flavor /5 in a beam that was born with flavor a 
and has travelled a distance L. This probability is given by 

P^p{L) = \A^p{L)\\ 

A^p{L) = 5]r/^»A,(L), (15) 

i=l 

where is the energy of the beam, and 



•4.W = i:5")=exp(.5-^), (16) 



Similarly, the probability for z/^ to oscillate into sterile neutrinos is 

3 oo 

PUL) = Y.Y.\Bam\\ (17) 

i=l k=l 



where 

oo 

5m(fe) = rE^°%'"exp 

n=0 

Of course, ELi Pc.p{L) + P„,(L) = 1. 



3 LED Sterile Neutrinos and SNO: Motivation and 
Choice of Parameters 

The flux of solar neutrinos has been measured by several collaborations using inde- 
pendent experimental techniques [jl4|. The observed flux is consistently lower than 



the predictions of the Standard Solar Model (SSM) by a factor of two to three, de- 
pending on the spectral sensitivity of the experiment. The most elegant explanation 
of this deficit is provided by the hypothesis of neutrino oscillations. The neutrinos 
produced in the Sun are in the flavor (z/e) eigenstate. If this state is a nontrivial su- 
perposition of the mass eigenstates, some of the neutrinos would change flavor on the 
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way to the Earth. Because the terrestrial experiments are predominantly sensitive to 
electron neutrinos, this effect leads to a suppression of the measured fluxes. 

Until recently, the experimental data could be explained either by an oscillation 
of z/g into other active species, z/^ or z/^-, or by its oscillation into a sterile neutrino 
sector, denoted by Ug. However, the best fit was for oscillations into active neu- 
trinos, especially when the Super-Kamiokande day-night and spectrum information 
was included 0. The case has recently been greatly strengthened. The SNO col- 
laboration has unambiguously demonstrated the presence of non-electron active 
neutrinos in the solar flux. Moreover, the total flux of solar neutrinos implied by SNO 
and Super-Kamiokande data is in excellent agreement with the SSM calculations |T5 |. 
These results mean that only a fraction (if any) of the solar neutrinos get converted 
into sterile species. 

In our framework, there are three active species, namely the left-handed brane 
neutrinos z/£. These states can mix both among themselves and with the excited KK 
states of the bulk neutrinos, z/^*-"'''. The bulk neutrinos do not carry any SM charges, 
and therefore are sterile. Thus, the SNO data indicates that the solar neutrino deficit 
is primarily due to the oscillations amongst the active brane states, and allows one 
to put constraints on their mixings with the sterile KK states. 

In the model proposed by Dvali and Smirnov [Q, the solar neutrino deficit is 
explained by the oscillations of z/£ into the KK excitations of the bulk neutrino of the 
same flavor. The other active flavors, z/j^ and z/J^, play no role. (In the language of 
section |2.2| , this corresponds to l^^ = 1, l'^'^ = l^^ = 0.) It is immediately obvious that 
this model is in contradiction with the SNO data: it predicts that no non-electron 
active neutrinos should be observed in the solar flux. To accommodate the data, it 
is necessary to introduce non-zero mixing among the active neutrinos. 

Motivated by the SNO data, we will take the following approach. We assume that 
in the limit in which the sterile KK neutrinos are decoupled {R limit) the oscil- 
lations amongst the remaining active species provide a good fit to the solar neutrino 
data. To be specific, we assume that the large mixing angle MSW solution to the 
solar neutrino problem, which provides the best fit to the data , is correct. We fix 
Sml^i = 6m\^ = (m^)^ — (mf )^ = 3.7 x 10~^ eV^ and take the mixing matrix elements 
to be Z"^^ = cos (3, = sin j3, with tan^ j3 = 0.37. (Varying these parameters over their 
allowed ranges does not significantly affect our conclusions.) Furthermore, we will as- 
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sume that the atmospheric neutrino anomaly observed by Super-Kamiokande [|T^] is 
also due predominantly to the oscillations between the active species (z/^ ^ z/,-). This 
assumption is supported by the analysis of the neutral-current-enhanced events, the 
zenith angle distribution, and r-appearance candidates in Super-Kamiokande |T8|| , 
and by the zenith angle distribution observed by MACRO [1^, which strongly dis- 



favor the possibility that the anomaly is due to ^ Vg oscillations. A good fit to 
the Super-Kamiokande data is obtained if we choose the mixing between z/^ and 
to be maximal, and the mass splitting l^m^j^l = |5m32| = 3 x 10^'^ eV^. Finally, we 
(generally) choose l"^^ to be because of the strong constraint on this matrix element 
provided by a combination of atmospheric neutrino data and reactor experiments. 
With these choices, the only remaining free parameters are the radius of the extra di- 
mension R and the absolute scale of the Dirac masses. For the latter, we will consider 
three possibilities: (1) the "hierarchical" mass scheme, 0, ^ m^, implying 

mf ^ (5"^atm)^ ~ 0.055 eV and mf ^ (<5mD^ ~ 0.006 eV; (2) the "inverted" 
mass scheme, ^ 0, mf ^ ^ (^"^atm)^ ~ 0.055 eV; and (3) the "degenerate" 
scheme^, with mf ^ ^ ^ 1 eV, where the common value has been chosen 
to be close to the cosmological limit of ~ 4.4 eV on the sum of the light neutrino 
masses [^. In the next section, we will present experimental constraints on the 
radius R in each of these three schemes. 

An important advantage of our perturbative approach is that it allows us to 
obtain simple analytic results throughout the analysis, making the underlying physical 
picture quite transparent. A more complete approach would be to fit the data allowing 
all the parameters of the model (3 Dirac masses, 3 mixing angles, 1 CKM phase and 
the radius of the extra dimension R) to vary, and then obtain bounds on R by 
marginalizing over the other variables. The resulting bounds would necessarily be 
less restrictive than the results of our present approach. However, as we explained 
above, current experimental data disfavors the possibility that oscillations into sterile 
states play a dominant role for either solar or atmospheric neutrinos. This implies 
that the results presented here should provide a good first approximation which could 
then be refined by a more systematic analysis. 



■''The original motivation for the degenerate scheme, i.e., mixed hot-cold dark matter models, has 
now largely disappeared, but it is still a logical possibility. Note also that this scheme requires a 
fundamental scale Mp larger than 100 TeV. 
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4 Experimental Constraints on R from Neutrino 
Oscillation Data 



Having stated our approach, we will use the formalism developed in the previous two 
sections to derive constraints on the radius of the extra dimension R from a variety 
of sources, such as reactor and accelerator searches for neutrino oscillations and solar 
and atmospheric neutrino flux measurements. We will next consider each of these 
constraints in turn. 



4.1 Reactor and Accelerator Experiments 

Large fluxes of anti-electron neutrinos are produced at nuclear power reactors. If the 
flux can be either predicted accurately or measured by a nearby detector, measuring 
the Ue flux at a certain distance L from the reactor gives the electron neutrino survival 
probability Pee{L). Several experiments, such as CHOOZ |^ and Bugey have 
utilized this approach to search for neutrino oscillations. Their results are consistent 
with no oscillation hypothesis, that is, Pee{L) — 1. 

The electron neutrino survival probability is determined by Eq. (pISl): 



Pee{L) 



2 



(19) 



j:\n'A{L) 

i=l 

where Ai is given by Eq. (0). For /^^ ~ and the values chosen for the mf, 
Pee{L) ^ 1 — Pes{L). FoT the case of small ^ that we will be primarily interested in, 

PUL) ^E^m' t( ^~ + oia (20) 



k=l 



where the oscillation phase is 

, ^ (A. -A )L ^ k_ ( D)2\ _L (21) 

It is convenient to evaluate the probabilities numerically, keeping only a finite number 
of terms in the sums in Eq. (|T^) or Eq. (PP]). (The oscillation probability into the 
/c-th KK state is suppressed by k'"^ for large k.) 

For the CHOOZ experiment, the distance from the reactor to the detector is 
L ^ 1 km, and the neutrino energies Ey range from about (1—6) MeV. The ratio of 
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the measured flux to the no-oscillation Monte Carlo prediction reported by CHOOZ, 
7^ = 1.01 ± 0.028(stat.)±0.027(syst.), implies^ 

Pee{L) > 0.942 [Pe. < 0.058] (22) 

at 90% c.l. Taking an average value L/E^ ~ 300 m/MeV, we obtain the following 90% 
c.l. constraints on the radius of the extra dimension: (1) in the hierarchical scheme, 
^2 < 0.43, corresponding to l/i? > 0.02 eV; (2) in the inverted scheme, ~ ^ < 0.13, 
1/R > 0.60 eV; and (3) in the degenerate scheme, 6 ~ 6 ~ ^3 < 0.13, 1/R > 10.9 
eV. We present all the corresponding upper bounds on R in Table |l|. The constraints 
in the inverted and degenerate schemes are quite a bit stronger than in the, perhaps 
more natural, hierarchical scheme. There are two reasons for this. Firstly, since 
^ = \/2mR, the same bound on ^ results in stronger bounds on R if the active 
neutrinos are heavy. Secondly, with the parameters of the CHOOZ experiment, the 
phases of the oscillations into the flrst few KK states are small in the hierarchical 
case, suppressing the oscillations. The higher modes have larger oscillation phases 
but small mixings. 

These results were obtained numerically using ([TBI) , but it is instructive to dis- 
cuss the approximate formula (pOf), which yields almost identical results, even in the 
hierarchical mass scheme. In the hierarchical case, 

PUL) - 2 sin^ P il £ (^^1?^) ■ (23) 

The phases 02 of ^^^t few KK states are not large, so the cos 02 terms must be 
kept. We have verifled that integrating over the neutrino energy Ey, which smooths 
the oscillations, gives results almost identical to using the average LjEy. In the 
inverted and degenerate cases, the 0j^2 large and the cos0j^2 average to zero. 
Since ~ ,^2 in these cases, one flnds Pes{L) ^ ,^|7r^/3. 

Our assumption = yields the most conservative limit. Relaxing this as- 
sumption would result in more stringent limits because 1 — Pee{L) would receive new 

§ These results and similar ones throughout the paper are obtained by the conservative method 
of assuming a Gaussian distribution for Pee{L) in the physical region < Pee{L) < 1, using the 
central value and uncertainty from TZ but renormalizing so that the total probability for finding any 
value of Pee{L) in the physical region is unity. 
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positive contributions ~ (2 + ^^3), where the two terms are respectively due 
to oscillations into active neutrinos and into the third KK tower. 

The analysis of the Bugey experiment is very similar. The results are presented 
in Table |l|. Because of the small distance between the reactor and the detector in this 
experiment, it is not possible to put significant bounds on R in the hierarchic scheme: 
oscillation phases are substantial only for very high KK modes, but the mixings of 
those modes with the active neutrino are severely suppressed. In the other two mass 
schemes, the bounds are similar to CHOOZ, although somewhat weaker due to a 
lower precision of the flux measurement. 

Neutrino oscillations can also be searched for at accelerators, where beams of 



muon neutrinos can be produced. The CDHS experiment |^ at CERN has searched 
for the disappearance of muon neutrinos from a beam with average energy ^ 3 
GeV, by measuring the ratios of the numbers of events in two detectors at 130 and 
885 m from the source. The negative results of this experiment can, in principle, yield 
bounds on R. However, for the L/E^ range of the experiment we obtain significant 
bounds only for the degenerate scheme. Averaging over the neutrino energies in the 
CDHS range and comparing with the measured ratio, we obtain the constraint that 
1/R must be greater than around 6 eV, with a small allowed window around 4 eV. 
Because of the imprecise nature of the degenerate solution, we just quote 1/i? > 4 
eV. 

More recently, the CHORUS and NOMAD ^ collaborations have searched 
for an appearance of tau neutrinos in a beam originally composed of v^. The non- 
observation of tau neutrino events has led to strong bounds on the probability of 

—>■ oscillations, of order P < 10^^. Unfortunately, this constraint does not put 
significant bounds on R. With the parameters of these experiments, the phases of the 
oscillations into the first 10-100 modes are small in hierarchical and inverted schemes. 
For higher modes, the oscillation amplitudes are of order ^ jr? ^ since the KK mode 
has then to oscillate back into to be detected. Thus, the corresponding probability 
is severely suppressed: P^r ~ ^'^/'^'^- The resulting constraints are of the order ^ ~ 1 
and are not competitive with bounds from neutrino disappearance experiments. For 
the degenerate scheme the oscillation phases could be substantial even for low-lying 
modes, but because of unitarity cancellations the resulting oscillation probabilities 
are severely suppressed, and no useful bound can be obtained. (The cancellations 
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appearing in this case are analogous to those discussed in section |^.) 

Another accelerator experiment, LSND, has reported evidence for neutrino oscil- 
lations. We postpone the discussion of this result until section |[ 



4.2 Atmospheric Neutrinos. 



The Super-Kamiokande experiment has measured the fluxes of electron and muon 
neutrinos produced in the atmosphere by cosmic rays. While the flux of electron neu- 
trinos is well described by a Monte Carlo (MC) simulation based on the no-oscillation 
hypothesis, the muon neutrino flux deviates significantly from this prediction. In 
particular, the flux of up-going muon neutrinos is suppressed in comparison with the 
down-going flux. The data, including the corresponding zenith angle distribution, 
is consistent with the hypothesis that on their way through the Earth the up-going 
muon neutrinos undergo oscillations that convert them into i/^'s. (An alternative pos- 
sibility of Ufj^ Ug oscillations is strongly disfavored |18[-) Fitting the data within the 
two-flavor scheme yields the mass difference |5^32l ~ 3 x 10~^ eV^, and the mixing 
angle ^23 ~ 7r/4 (maximal mixing.) 

For the case we are considering, the oscillation of i/^^e into sterile KK states will 
involve large oscillation phases and therefore lead to a suppression of the overall 
Uf, and/or z/^ fluxes that is independent of the zenith angle. The Super-Kamiokande 
collaboration have included the possibility of such L/£',^-independent flux changes 



in their analysis to take into account the theory uncertainties in the absolute and 
relative flux calculations. However, their results can also be used to constrain L/ Ey- 
independent oscillations into sterile states, provided the theoretical uncertainties in 
the fluxes are properly included. 



The best-fit values of the Super-Kamiokande data [Tj] , assuming a, anal- 
ysis, suggest that there is a slight overall excess in the ratio of the measured to the 
central value of the simulated number of events. The data also yield a deficit for 
the ratio of the /x-like to e-like events, assuming the same two-flavor analysis. In 



the analysis of the data |T^, the overall excess and the deficit mentioned above are 
parameterized by a = 0.034 ± 0.25 and (3s = —0.059 ± 0.08, respectively, where the 
observed sub-GeV z/g and rates relative to the number expected for no oscillations 
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into sterile states are 

w = (i + «)(i + y)- (24) 

The subscript s refers to the fact that this value is obtained from the sub-GeV events, 
which yield the most stringent constraints. We have assumed that the errors are 
dominated by the theoretical flux uncertainties. Quantitatively, the experimental 
results imply 

- a = i[P(z/^ z/,) + P(z/e ^ i^s)] (25) 

and 

= P(z/e^z/.)-PK^z/,). (26) 

In what follows, we use the one-sided 90% c.l. upper bounds on —a and ±Ps to derive 
our constraints. 

In the hierarchical mass scheme, the active neutrinos predominantly mix with the 
z/^*^^\ since ^3 ^ ^2? ^1 ~ 0. Since /^'^ = 0, this mixing can only occur for the muon 
neutrinos, decreasing their flux compared to the best-fit value; the electron neutrino 
flux remains unchanged. Eq. (p6D implies that at 90% c.l., 

Pii^^ i^s) - P{i^e ^ i^s) < 0.17. (27) 

Thus, ignoring the z/g oscillations, and in the limit 03 ^ 1, we obtain P^u^ i^s) = 
^IvrVe < 0.17, which yields 1/P > 0.24 eV. 

In the inverted mass scheme with ^1 ~ ^2 = ^ and ^3 ~ 0, two KK towers, u]^^''^ 
and contribute to the oscillations. Both electron and muon neutrinos can mix 

with these states, so both fluxes are reduced. It turns out that the suppression is 
more significant for the electron neutrino, so we derive our constraint from the upper 
bound on (3s- Again, assuming the phase energy averaging, we find P(z/e t^s) = 
2P(z/^ Us) = ^^2/3. Thus, ^2vr2/6 < 0.1, which yields 1/P > 0.32 eV. 

Finally, in the degenerate mass scheme, with ^i~^2~^3 = ^;all three KK towers 
can mix with the active neutrinos with an approximately equal strength. In this case, 
the oscillations into the KK states suppress the electron and muon neutrino fluxes by 
the same amount. In this case, we use Eq. (^) to flnd (1/2) [P(z/^ —>■ Us) + P{ue 
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Us)] < 0.39, at 90% c.l. We have P{u^ u^) = P{ve ^ v,) = where phases 

are energy averaged. This yields .^^vr^/S < 0.39 and hence 1/R > 4.1 eV. As we 
can see, the bounds from the atmospheric data are quite stringent. These results are 
summarized in Table 0. 



4.3 Solar Neutrinos 

As we already discussed in section H, the results of the SNO collaboration P make it 
unlikely that oscillations into sterile KK neutrinos play a major role in explaining the 
solar neutrino deficit. In particular, combining the SNO and Super-Kamiokande |^ 
results, and using the new SSM estimate of the flux of neutrinos [|l^] that incor- 
porates a recent more accurate measurement of the '^Be{p,'jYB reaction [16], one 
obtains 

J2 P{^e l^a) = 0.62±0.21, P(z/e ^ z/,) = 0.083±0.21, P(z/e ^ l^e) = 0.30±0.05 

a=fi,T 

(28) 

for the probabilities for an initial ^B neutrino in the SNO/Super-Kamiokande energy 
range to oscillate into z/^, z/^, or a sterile state, or to remain a Ue, respectively. This 
leads to 

Pii^e i^s) < 0.40 (29) 

at 90% c.l. (The corresponding limit would be 0.41 using the older cross sections.) 

To compute the fluxes of solar neutrinos, we need to take into account matter ef- 
fects in the Sun. In vacuum, the flavor-basis active neutrino states can be decomposed 
into the mass eigenstates according to 

3 oo 
i=l n=0 

(This is just Eq. (|l^) of section ^12] ; we have dropped the subscript L to avoid clut- 
tering and defined uj^^ = /"*L°".) In matter, an additional, flavor-dependent effective 
mass term is generated for active neutrino species. The decomposition ( |5D| ) is modified 
by this effect: 

3 oo 

^a = EEf^S^"^^^"^ (31) 

i=l n=0 

where z/J^"-* and ?7™ are respectively the eigenstates and mixing matrix elements in 
matter. These of course depend on the matter density, and therefore change as the 
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neutrino travels through the Sun. We assume that neutrino propagation is adiabatic, 
that is, the fraction of each of the mass eigenstates in the beam remains constant. 
This is an excellent approximation for the large mixing angle (LMA) solar neutrino 



parameters that we are using [14, 26|. Since all the phases associated with solar 



neutrino oscillations are large, we can neglect the interference effects between the 



matter eigenstates. Then, the oscillation probability is given by a simple formula 



i,n i,n 

In this equation and below, i/^"^ and f/^*'"^ are evaluated at the neutrino production 
point in the core of the Sun, and we have assumed that the neutrino flux is measured 
in the vacuum (as is the case for solar neutrinos.) 

The values of the mixing angles f/™*'"^ in (|3^) can be computed exactly if the 
neutrino energy and electron and neutron densities in the region of the Sun where the 
neutrinos are created are known. This in turn depends on the number and properties 
of the Mikheyev-Smirnov-Wolfenstein (MSW) resonances passed by the neutrino 
on its way to the vacuum. 

For the energies of the neutrinos observed by Super-Kamiokande and SNO, an 
MSW resonance can only occur for mass-squared splittings less than about 10~^ eV^. 
Hence, in our case, resonant conversion is only possible between the two zero-mode 
mass eigenstates, u^^^^ and z>^*^°\ The higher mass states, which are mainly the sterile 
KK excitations, do not undergo resonant conversion and are little affected by matter 
effects for the masses required by the reactor and atmospheric data. Thus, there 
is only one resonance to consider. We will assume that all solar neutrinos observed 
by SNO and Super-Kamiokande cross this resonance, and that for all of them the 
propagation through the resonance is adiabatic, which are excellent approximations 
for the LMA parameter region. The electron neutrino mixing angles in the core of 
the Sun are then 

= 0, = If/iJ^I^ + \ui?\\ (33) 

with all the other matrix elements f/™*^"^ being equal to their vacuum values. 

Formulas (|32|) and (p3|), together with (0) and the small-^ expansions of the 
vacuum mixing angles (ll^), allow us to estimate the sterile component of the solar 
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Experimental Bounds 


Experiment 


Hierarchical 


Inverted 


Degenerate 




(cm, eV) 


(cm, eV) 


(cm, eV) 


CHOOZ 


(9.9 X 10-^0.02) 


(3.3 X 10-5,0.60) 


(1.8 X 10-6,10.9) 


BUGEY 


none 


( 4.3 X 10-^0.46) 


(2.4 X 10-6,8.3) 


CDHS 


none 


none 


(5 X 10-6,4) 


Atmospheric 


(8.2 X 10-^0.24) 


(6.2 X 10-^0.32) 


(4.8 X 10-6,4.1) 


Solar 


(1.0 X 10-3,0.02) 


(8.9 X 10-^0.22) 


(4.9 X 10-6,4.1) 



Table 1: Upper bounds on R (cm) at 90% c.l. and the corresponding lower bounds 
on 1/R (eV) from various measurements. 

neutrino flux. For example, in the hierarchical mass scheme, 

Piiy,^us) = l- E P(^^e^^^a) = Y(l + sin'/5)e2- (34) 

The experimental bound P{ve — > z/^) < 0.40 at 90% c.l. then implies 1/R > 0.02 
eV. The analysis in the other two mass schemes is analogous except that P(z^e 
Us) = 7r^^|/3, implying R > 0.22 eV (inverted scheme) and 1/R> 4.1 eV (degenerate 
scheme). These results are shown in Table ^ 



5 The LSND Experiment 

The LSND experiment has searched for 17^ —>■ z7g oscillations and reports an excess of 
events, corresponding to an oscillation probability P{T'f_i z/g) = (0.264 ± 0.067 ± 
0.045)% |TD[. Similar results are obtained for z/^ u^.. The data suggests that the 
oscillations occur for 0.2 eV^ < < 10 eV^, with the lower Sm^ values favored by 



the nonobservation of an oscillation signal in the KARMEN experiment p^. Thus, 
it seems that mass scales of ~ 1 eV must be present in any model that attempts to 
describe the LSND results. In particular, one must introduce a fourth neutrino to 
account for LSND as well as the solar and atmospheric results, and this must be a 
sterile Ug due to the constraint on the number of light active neutrinos from the Z 



width p8l. Many authors have analyzed such four-neutrino schemes [29|. The best fits 
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are obtained in the so-called 2 + 2 schemes, in which there are two closely spaced pairs 
or states, with the pairs separated by ~ 1 eV and a small mixing between the pairs 
to account for the LSND results. However, the simplest versions, in which the pairs 
consist of either (a) {ue, ^s) and (z/^, Ur), or (b) {u^., Vr) and (z/^, Vg)^ are now excluded 
by the solar and atmospheric data, so that the partners of the and z/^ would have 
to be mixtures of Vg and Vr- The alternative 3 + 1 (or 3 + p, with p > 1) scheme 
involves three closely-spaced active neutrinos (similar to the hierarchical, inverted or 
degenerate schemes), separated from one or more mainly-sterile states by about 1 eV. 
i^pL oscillations are therefore a sub-leading effect involving the mixing of both 

and z/g with the sterile state. Such schemes are strongly disfavored because it is 
difficult to obtain a large enough effect for LSND while still respecting the reactor 
and accelerator z/g and disappearance constraints, but are still barely possible. 

Although the (3, 3) model contains KK states of mass ~ 1 eV in the hierarchical 
scheme for large mode numbers k ~ 100, transitions among the active states are 
suppressed by the smallness of (so far set to zero) and bvn? for the active states, 
as well as by the 1/fc factors in the mixing. Even allowing |/^^| «i 0.2, calculations 
based on our formalism still fall an order of magnitude short of the central value for 
-P(j^Ai ^e) quoted above. The inverted and the degenerate schemes result in roughly 
the same or even smaller probabilities. In generic 3 + 1 or 3 + p schemes the rate is 
simply proportional to the fourth power of mixing angles Q , i.e., P{y^ — > Vej oc O^^O^^. 
However, in our model there are cancellations between the two KK towers so that 
P{V, ^ Z7,) oc lel - e?P = 4|<5mLi?T (or m6ml,^R'\^ for l^^ ^ 0). Thus, 
there are no enhancements from the larger Dirac masses in these schemes, and also 
R cannot be larger than in the hierarchical case because of other constraints. 

Hence, it seems that an extension of the (3, 3) model is needed in order to address 
the LSND results. Several authors have considered the introduction of Majorana 
masses of unknown origin for the brane states. In the following, we propose two 
other extensions that seem to provide the necessary ingredients to accommodate the 
LSND data. We give only simple estimates of the effects. To determine whether these 
extensions can indeed provide a consistent explanation of all neutrino oscillation data 
requires a more careful study, which is outside the scope of the present work. 
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5.1 The (4, 4) Model 

Consider two extra sterile neutrinos, denoted by A^* and u'l, residing in the bulk and 
on the brane, respectively; z/£ is left-handed. A mass term that couples these two 
states is consistent with all the symmetries of the theory. It has the form 

C^ass = XMFTyliy%{x,0) + h.c., (35) 

where A is a coupling analogous to Xap in Eq. (0), and z/|j is the right-handed compo- 
nent of N^. Note that Mp is the only natural choice for the mass scale. Introducing 
a dimensionless Yukawa coupling h = XMp and performing the KK decomposition 
of z/^, we obtain mass terms of the form 



[vTui + V2 E 17? Vij + h.c, (36) 
where the mass scale is given by 

ms = h^. (37) 

M PI 

For Mp ~ 100 TeV and /i ~ 1, this mass scale is in the 1 to 10 eV range. Thus, the 
scale required to account for the LSND results has been obtained without any fine- 
tuning, using only the natural parameters of the theory. Note that the symmetries 
of the theory also allow the brane Yukawa couplings between z/^ and the active left- 
handed neutrinos z/£, which will lead to mixings between the active species and the 
fourth sterile KK tower. It seems plausible that the solar, atmospheric and LSND 
results could be accommodated simultaneously in this framework; however, a more 
careful analysis is clearly called for to verify this. 



5.2 Two Large Dimensions 

As discussed above, the small rate expected for LSND in the canonical (3,3) model is 
due to a cancellation between the contributions of two KK towers, so that even in the 
inverted and degenerate schemes the rate is proportional to |^2~^iP l^'^^Pl^l P 
rather than |^2|^ or j/'^^Pl^sl'^. This cancellation can be traced to the fact that all of 
the bulk neutrinos propagate in the same large dimension. The cancellation could be 
broken in the case in which there are two or more extra dimensions of unequal radii. 
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(3^ 



In the following, we consider a simple phenomenological model in which the three 
bulk neutrinos each propagate in only one extra dimension, of radii -Ri, R2, and i?3, 
respectively, where two or more of the Ri may be different. We make no attempt to 
construct a full model or derive all of the constraints, for our goal is only to illustrate 
a possibility. 

One can then show that in the limit when the 1/ Ri are all large compared to the 
Dirac masses mf , the amplitude A^^i^L) in (|T3]) becomes 

J fie 

where S = A/2m^R, is the analog of the Dirac mass matrix in (^, R is the 3x3 
matrix diag i?2, -R3), and the matrix P is given by 

k=i 

with (f)^ = [{k'^L) / (2i?;y)]R~^. Writing = l^^mj^r, where = diag (mf , m^, m^) 
is diagonal, we see that Af^^iL) does not vanish for ^ provided that the left 
and right unitary matrices / and r are not equal. As a simple two-family example 
with / and r given by rotations with angles Pl and Pr, respectively, one finds 



Pn--jRi-[P..--jRl 



Af,e{L) = 2(mf )^ {-slCr + cl-sr) {clCr + sl-sr) 

L \ ^ / \ ^ / 

(40) 

where sl,r = sin/J^jj, and similarly for cl,r- This is easily generalized to the three- 
family case that can accommodate atmospheric neutrino results. Thus, it is clear 
that the z/g ^ oscillation rates can be enhanced compared to the (3, 3) model. 
However, a full investigation of this proposal is beyond the scope of this paper. 



6 Conclusions 

In this paper, we have studied neutrino oscillations in the context of theories with 
LED's. In particular, we focused on a 5-dimensional model with three active brane 
and three sterile bulk neutrinos, coupled to the 4-dimensional Higgs field. This setup 
yields small Dirac masses for the active species, in addition to a tower of sterile KK 
neutrinos. An attractive feature of this (3, 3) model is that it does not introduce 
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any ad hoc parameters, such as small 4-dimensional Majorana masses. With an extra 
dimension of radius R, the lightest KK neutrinos have masses ~ Our approach 
was based on the assumption that a three-flavor analysis of the data provides the 
best fit for the observed oscillations, and we have treated the effect of the LED as 
a perturbation. Hence, in our treatment, the Dirac masses are taken to be much 
smaller than 1/R. 

The recent solar neutrino data from the SNO collaboration, in conjunction with 
the results of the Super-Kamiokande experiment for atmospheric neutrinos, signifi- 
cantly constrain the size of the effect of sterile states on the solar neutrino deficit. In 
particular, models based solely on active-sterile oscillations, as in Ref. are now 
ruled out. This fact motivated us to determine what bounds the new data place on the 
contributions of KK sterile states to neutrino oscillations. Since these sterile states 
originate in the extra dimensions, we obtained constraints on the size R of the largest 
of the extra dimensions. In addition to the new solar data analysis, we also performed 
analyses of reactor, accelerator, and atmospheric data. The (3, 3) model studied here 
cannot accommodate the LSND data. However, we proposed two natural extensions 
of this model to address the LSND results. One is a (4, 4) model with two extra 
sterile neutrinos, one on the brane and another in the bulk, which can yield a 3 + 1, 
or, less naturally, a 2 + 2 four-neutrino scheme. These are phenomenologically similar 
to four neutrino schemes obtained in other theoretical frameworks, yielding equally 
good or poor descriptions of the data. However, there are no Majorana masses, so 
there is a conserved lepton number and neutrinoless double beta decay is forbidden. 
The other extension uses two LED's of unequal radii, in each of which a bulk neutrino 
is sequestered. We provided estimates of the effects in these two cases and concluded 
that the LSND results seem to be accommodated in our proposals, although a more 
careful analysis is warranted. 

In this paper, we have concentrated on constraints from laboratory, solar and 
atmospheric neutrino oscillation experiments. There are also complementary con- 
straints from cosmology and astrophysics. In particular, the degenerate scheme may 
be excluded by considerations of big bang nucleosynthesis, cosmic microwave back- 
ground radiation, and diffuse extra-galactic background radiation [^, but the hier- 
archical and inverted schemes are consistent. There are also potentially very strong 
limits on all of the schemes coming from the constraints on the energy carried away by 
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the sterile KK modes from deep inside a supernova core. Qualitative estimates 
concentrating on a single family suggest that the constraints may be considerably 
stronger than those in Table 1. A more detailed study would be very useful. Bounds 
on the model considered here can also be derived from studying the effects of KK 
neutrinos in electroweak processes 0. 

Our main results are contained in Table |1|. These are to be interpreted as con- 
straints on the size of the largest of the extra dimensions, regardless of their total 
number. If, as seems plausible, the pattern of masses in the neutrino sector is like 
those of the other fermionic sectors, and hence hierarchical, the strongest bound comes 
from the atmospheric neutrino data and requires R to be less than about 0.8 fim. The 
bounds in the inverted and degenerate mass schemes are even stronger. On the other 
hand, the currently discussed Cavendish type experiments have sensitivities only of 
order 50 /im, and an improvement of two orders of magnitude would be necessary to 
reach the region allowed in our model. Thus, to the degree that the (3, 3) model is a 
natural context for small Dirac masses in theories with LED's, our analysis suggests 
that neutrino oscillations seem to rule out the possibility of observing the gravitational 
effects of the extra dimensions in the foreseeable future. 
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